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ABSTRACT 

I consider a toy model of self-regulated black hole accretion. The black hole grows 
through Bondi accretion and a fraction of the accretion power is distributed as thermal 
feedback into the surrounding gas. The gas expands or contracts until AGN heating 
and radiative cooling balance each other. The balance of heating and cooling is used 
to determine a quasi-equilibrium temperature at which the black hole accretes in self- 
regulated equilibrium with the surrounding intergalactic medium. This temperature 
grows with the black hole mass. The temperature increase is very steep around a 
critical black hole mass due to the shape of the cooling function. The quasi-equilibrium 
temperature cannot exceed the virial temperature or the AGN will drive a thermal 
wind. This limits the black hole mass to a maximum value determined by the depth of 
the potential well. In the regime in which cooling is dominated by bremsstrahlung, this 
model determines a relation between black hole mass and halo characteristic velocity 
of the form M, oc v A . The predictions of the model are consistent with the observed 
black hole mass - bulge velocity dispersion relation. 
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1 INTRODUCTION 

All massive spheroids for which the data are sufficiently 
good show evidence for a central black h ole. The rel ation 
between black hole mass and bulge mass iMarconi fc Hunt] 
120031: lHaring fe RbdEool) suggests that the black hole con- 
tains ~ 0.2% of the bulge stellar mass. The gravitational 
binding energy released by matter accreted onto the black 
hole is ~ 10% of the rest mass energy. A small fraction of this 
energy is sufficient to affect the dynamics and the thermal 
state of the intergalactic medium (IGM) dramatically. 

Active galactic nucleus (AGN) heating was in- 
troduced t o solve the cool ing flow pr oblem in galaxy 
clust er s (iTabor fc Binnevl Il995t iBinnev fc Tabol 
[ 19951: iTucker fc Davidl 119971 ICiotti fc Ostrikerl Jl997t 
iGwalier^^tHnl l2002t te^zkowski^^B^gelmanT I2OO2T1 
and the related problem of the e n tropy floor 
jRovchowdhurv et alJ I200I lLapi et aljl2005l) . ISilk fc Reesl 
(1998) went further and proposed that AGN outflows are 
an essential part of the galaxy formation process. In this 
scenario, it is AGN feedback that determines the observed 
relation between black hole mass and velocity dispersion of 



and lDi Matteo et alJ ]2005h have developed more elab orate 
versions of the original idea by ISilk fc Reesl (^,998) that 
AGN winds terminate star formation in the host galaxy 
(the reference list does not pretend to be complete). 
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The interaction of the black hole with the IGM is a 
com plex problem, which requires hydrodynamic simulations 
fe.g. | Quili s^^lj20pj]:[Revnold^et all200jJ:[CJrurazov^^^J 
" 20021 hlevnolds et all 120021: iBasson fc Alexander! 
2003: lOmma fc Binnevl l2004t iRuszkowski et al] 12004 
Briiggen et al] l200Et iDi Matteo et al] 120051: IZanni et al] 
2005), but here I want to concentrate on a simple aspect, 
where an analytic approach is possible and can offer a 
better insight into the results of hydrodynamic simulations. 

I consider a black hole at the centre of a profile with 
a core, where the gas is approximately homogeneous and 
isothermal. The black hole accrete mass and returns heat to 
its surrounding environment. The core expands or contracts 
until AGN heating and radiative cooling balance each other. 
The balance of heating and cooling determines the relation 
between the mass of the black hole and the temperature of 
the IGM at the equilibrium. If the equilibrium temperature 
is much larger than the virial temperature, the system will 
only admit outflow solutions. 
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2 MODEL OF THE BH - IGM COUPLING 

A black hole of mass M. surrounded by gas with density p 
and speed of sound c s accretes at a rate determined by the 
iBondi Jl952h formula: 



M. = 4Tva(GM.) 2 4 
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where G is the gravitational constant and a is a number of 
order unity. In the second equality, the speed of sound c s is 
rewritten in terms of the temperature T. Here 7 is the adia- 
batic index, k is the Boltzmann constant, m p is the proton 
mass and fj, is the mean particle mass as a fraction of the 
proton mass. Eq. Q assumes that the accretion is spherical 
and that the gas is approximately static at infinity, but is 
widely used to es timate M, even when t hese conditions are 
not satisfied fe.g. lDi Matteo et al.ll2005T) . 

A fraction j3 of the power generated by the accretion of 
matter onto the black hole heats the surrounding gas at a 
rate Qheat- If the heat is distributed over a gas mass M gas , 
then the heating rate per unit mass is: 
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where c is the speed of light and e ~ 0.1 is the energetic 
efficiency of black hole accretion. In the second line, where 
we insert the numerical values, eo.i = e/0.1, 75/3 = 37/5, 
jUi6/27 = 27/i/16, X3/4 = 4X p /3 (X p is the hydrogen baryon 
fraction), M. 8 = M./1O 8 M , M gas8 = M gas /1O 8 M and 
T 6 = T/10 6 K, while riH is the number density of hydrogen 
atoms (free protons for an ionised gas). 

The cooling rate per unit mass of the heated gas is: 
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where An^ is the radiated power per unit volume. A is a 
function of the temperature T and the metallicity Z, and 
A_22 = A/10 -22 erg s _1 cm 3 . In the next Section I proceed 
to examine under what conditions the heating term and the 
cooling term balance each other off. 



3 LOCAL M. - T EQUILIBRIUM RELATION 

Both the heating (Eq. |5J and the cooling (Eq. [3J| rate scale 
linearly with the density of the gas. The difference is in the 
temperature scaling. Fig. 1 compares the the rates of heat- 
ing and cooling as a function of temperature. The curves 
are computed fro m the collisional ionis a tion equilibrium 
cooling function of ISutherland fc Dopital l|l993f) and show 
Q coo \{T) / (nuM^s) for different metallicities. The diagonal 
straight lines show the heating term Qh oa t (T) / (nrr Af gas ) and 
are parameterised by the black hole - IGM coupling factor 
£ = a/WM^ g/Mgas 8 (Eq.|5J, where S is the ratio between the 
gas density at the Bondi radius (in Eq.|2J and the gas density 
at the radius that contains a gas mass M gas (in Eq.^. 

At very low T, the heating term always dominates and 
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Figure 1. Qh ca t(T, ()/(n H M gas ) (diagonal lines) and 
Qcool(T, Z)/(nn Mg as ) (curves) for different values of the 
black hole - IGM coupling constant f and the metallicity Z. 



the temperature increases. At very high T, the cooling term 
dominates and the temperature goes down. In both cases, T 
will converge to an equilibrium temperature T cq determined 
by two parameters: Z and f. 

For each metallicity, there is a value of £, which we 
call £(Z), where the equilbrium extends over a temperature 
range (Fig. 1). For zero metals, this value is C(Z) ~ 1CF 3 
and the temperature range is 10 5 K < T Gq < 3 x 10 5 K. For 
Z = 0.1 Zq, (,(Z) ~ 0.025 and the temperature range is 
2 x 10 5 K < T cq < 10 B K. For Solar abundances, £(Z) ~ 0.1 
and the temperature range is 2 x 10 5 K < T cq < 3 x 10 6 K. If 
the parameter £, which measures the effective efficiency of 
AGN heating, is £ < Q{Z), then cooling dominates at all but 
the lowest temperatures and AGN feedback is negligible. 

For £ 3> C(-^)j cooling is dominated by bremsstrahlung. 
In this regime, the cooling function can be approximated as 
A_22 ~ 0.2vT 6 , so that equating ||5J and © gives: 

M.g 
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The equilibrium temperature is proportional to the black 

hole mass. The constant of proportionality depends on the 

black hole accretion efficiency with respect to the Bondi rate, 

on the energetic efficiency of AGN heating and on the gas 

mass on which the feedback is distributed. 

The equilibrium described in this Section is a local 

quasi-equilibrium. It is a quasi-equilibrium because T eq is 

not constant with time. The black hole accretes at the rate 

specified by Eq. Q and T eq must constantly increase to 

keep up with the increasing black hole mass (£ oc M?). 

However, substituting Eq. Q into Eq. Q shows that the 

1/2 

quasi-equilibrium accretion rate is M, oc M, , while the 
accretion rate without self-regulation would be M, oc Ml 
until the Eddington limit is reached. This model predicts 
that in the self-regulated regime M,/MEddingtou oc M^ 1 ^ 2 . 
The equilibrium is local because we made no assumption for 
the scale on which the feedback is distributed, which may be 
as small as the nuclear region or as large as a galaxy cluster. 
In the next Section we shall see how relating the feedback 
scale to the properties of the host system can generate a 
limit mass for black hole growth. 
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4 THE Mm - Tvir EQUILIBRIUM RELATION 
IN GALAXY CLUSTERS 

If AGN heating is relevant not only for the self-regulation 
of the black hole and for the gas in the central ~ 1 kpc, but 
also to keep the cluster gas hot, then the feedback must be 
distributed on a scale of the order of the radius of the cluster 
core (~ 100 kpc). In this scenario, one assumes that T cq ~ 
T v i r and that M gas is equal to the mass of the gas in the 
core of the hot gas distribution, and one uses the equation 
Qheat(T)/(n H M gas ) = Qcooi(r)/(n H Afgas) to determine the 
equilibrium black hole mass. 

IKomatsu fc SeliakI i200ll) have assumed that the dark 
matter distribution i n a virialised halo is described by the 
INavarro et alJ 1)19971) profile and have used this profile to 
compute analytic hydrostatic equilibrium solutions for the 
baryonic component (I verified with computer simulations 
that these are truly hydrostatic solutions). They contain 
three parameters: the virial mass M v i r , the virial radius r v i r 
and the mass fraction in hot gas, which, for a cluster, is of the 
order of the cosmic baryon fraction (~ 0.1). The first two are 
not independent when a critical density contrast (a redshift 
of collapse) is specified. The relation between the mass of the 
black hole M. and the virial velocity v v [ r — GM v i r /r v ir is de- 
termined from Eq. J1J, wh ere T eQ ~ GMm D M v i r /(3fcr v i r ) and 
Mgas is determined by the lKomatsu fc SeliakI l)200l|) model. 

In an isothermal sphere the virial velocity « v ir and the 
velocity dispersion a are related by v v i T = a\/2. Fig. 2 shows 
the M, - a relation between the mass of the black and the 
velocity dispersion of the host galaxy if the stars in the bulge 
of the central galaxy have the same velocity dispersion as the 
dark matter. The two lines correspond to different redshifts 
of collapse of the dark matter halo. The redshift of collapse 
z c is used t o calculate the virial over density with the fitting 
formulae of iBrvan fe NormarJ il998f) and thus to compute 
r v ir as a function of M v i r . The data points show the observed 
black hole mass - bulge velocity dispersion relation. 

The model and the data overlap already in the simple 
case a ~ (3 ~ 8 ~ 1 (Fig. 2). One must exert great caution in 
relating the virial velocity of the cluster « v jr to the velocity 
dispersion of the central galaxy a. A more careful analysis 
would show that D v ir/\/2 is an overestimate of the velocity 
dispersion of the central galaxy, which is compensated by 
an underestimate of 5, since the gas density is higher at the 
centre due to the presence of the central galaxy. 

I will illustrate this point by considering the special 
case of the radio source in M87, at the centre of the Virgo 
cluster. If we assume that the Virgo cluster has a mass of 
Mvir ~ 1.5 x 10 14 Mq and has virialised at low redshift, 
then v v i T ~ 690kms _1 and T V1I ~ 10 7 K. If we also assume 
that the baryons are in the form o f a ho t gas described 
by the model of IKomatsu fe SeliakI i200ll) , then the pre- 
dicted values for the central electron density and the mass 
of the baryons in the cluster core are n c ~ 0.004 cm~ 3 
and M c ~ 1.8 x 10 12 Af Q (in reality they will both be 
lower because some of the baryons are in stars). The virial 
temperature T v i r is consistent with the X-ray tempera- 
ture in the core of the Virgo cluster (e.g. IDi Matteo et all 
2003), but « v ir/\/2 ~ 500 km s - 1 is much larger than the 
a ~ 375 kms -1 determination bv lTrema ine et al.l (120021) for 
the stellar velocity d i spersi on in M87. lHarms et al.l (Il99l) 
and IMacchetto et alJ Jl997t) have estimated that M87 con- 
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Figure 2. The predicted M» - v v i T /\/2 relation (lines) compared 
with the observed M» - a relation (points with error bars) , where 
a is the velocity dispersion of the host bulge. The model is for 
bremsstrahlung cooling (Eq. |1J. The upper line corresponds to 
a redshift of collapse of z c ~ 2. The lower line corresponds to 
z c ~ and thus to a lower halo density. The squares with erro r 
bars are the black hole mass estimates oflTremaine et ahl |200~' 



The d iamonds with error bars are those of Fcrrarcs e fc Merrittl 
(2000). 



tains a supermassive black hole of M. ~ 3 x 10 9 Mq . 
IDi Matteo et all i2003T) have used this estimate for the black 
hole mass together with their X-ray observations of M87 and 
have shown that the mechanical power of the jet in M87 can 
be understood if the black hole accretes the surrounding 
hot gas at the Bondi rate (a ~ 1) and most of the accretion 
power is released mechanically (/3 ~ 1). They have also ob- 
served that the central de nsity of the gas in M87 is higher 
than our estimate from the lKomatsu fc SeliakI i200lfl model 
by a factor of 5 ~ 30 — 40. With these values, Eq. Q gives 
Mgas ~ 2 x 10 12 Afo, which is indeed consistent with our 
estimate for the mass of the baryons in the core of the Virgo 
cluster. 

With this normalisation issue in mind, the comparison 
in Fig. 2 demonstrates that this model can reproduce the 
type of relation observed in the data. The predicted M, - a 
relation is not exactly a power law, but is consistent with a 
power law relation of the form M. oc a v with r\ ~ 4. 



5 DISCUSSION AND CONCLUSION 

I have considered a very simple model of self-regulated black 
hole accretion. The black hole grows through spherical Bondi 
accretion and a fraction of the accretion power is distributed 
as thermal feedback into the surrounding IGM. The heating 
rate per unit mass of the gas depends on three parameters: 
the accretion efficiency with respect to the Bondi rate, a, 
the fraction of the accretion power converted into heat, /3, 
and the mass M gas on which the heat is distributed. If Af gas 
is much larger than the mass within the Bondi radius and 
the interaction between the AGN and the IGM happens at 
large scales (e.g. ~ 100 kpc), then one must also deal with 
the complication that the gas affected by feedback has a 
different density from the gas that feeds the black hole. 

If the black hole is not producing enough heat or if the 
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heat is distributed on a very large scale, the thermal energy 
input is radiated almost immediately and the equilibrium 
temperature is not much higher than the IGM temperature 
without AGN feedback. For a given metallicity, there is a 
critical heating rate per unit mass at which the heating and 
the cooling rates are comparable over a broad temperature 
range. This is the range where the cooling function decreases 
with the temperature. Above this critical rate, the gas is only 
able to cool at high temperatures, in the bremsstrahlung 
regime. For a given set of parameters, the critical heating 
rate corresponds to a critical black hole mass. Feedback is 
inefficient below this critical black mass, while it can rapidly 
heat the gas to very large temperatures after the black hole 
has passed this threshold. 

In the Bondi model without self-regulation, the black 
hole grows with ~ M,. Self-regulation reduces the power 
with which the accretion rate depends on the black hole 
mass (~ M^' 5 for cooling by bremsstrahlung). 

If the equilibrium temperature is higher than the virial 
temperature, the black hole drives a thermal wind and the 
approximation that the heating is quasi-static breaks down. 
This blow-out condition implies a maximum M. (T v ir), which 
gives rise to a relation of the type M. oc u* ir . 

The question is the relevance of the approximations 
made by the model: 

i) Bondi spherical accretion in a homogeneous medium 
supported by thermal pressure may be a bad estimate of how 
Mm depends on the black hole environment, particularly if 
the gas that feeds the black hole is cold and clumpy. 

ii) In a multiphase IGM, the phase which feeds the black 
hole may not be the same that receives the heat. One can 
imagine a scenario in which the AGN is fuelled through the 
accretion of cold clouds, but most of the AGN heating goes 
to the dilute hot gas that fills the space between the clouds. 
This is the most likely picture in powerful AGN triggered 
by galaxy mergers. 

iii) Feedback may be mechanical and only thermalised 
on very large scales after escaping from the galactic core 
in a collimated outflow. In our model this corresponds to a 
large value of Af gas , but the difference between the accre- 
tion scale and the thermalisation scale may determine an 
oscillation pattern rath er than an equilibrium solution (e.g. 
lOmma fc Binnevll20o4) . 

iv) The onset of a thermally driven wind invalidates the 
assumption that heating is quasi-static. In the presence of 
outflow rather than static boundary conditions, the Bondi 
form ula greatly overestimat es the black hole accretion rate 
(e.g. iDi Matteo et al]|2003l) . This will almost certainly be 
the case if the quasi-equilibrium temperature is much higher 
than the virial temperature of the system. 

The toy model presented in this paper is mostly relevant 
to low power AGNs fed through the accretion of hot gas in 
the core of massive haloes, but it provides a general rule of 
thumb to determine under what conditions AGN feedback 
is important and self-regulation is possible. 



6 ACKNOWLEDGEMENTS 

I acknowledge useful conversation with Arman Khalatyan 
and Romain Teyssier. 



REFERENCES 

Basson J. F., Alexander P., 2003, MNRAS, 339, 353 
Binney J., Tabor G., 1995, MNRAS, 276, 663 
Bondi H., 1952, MNRAS, 112, 195 

Briiggen M., Ruszkowski M., Hallman E., 2005, ApJ, 630, 
740 

Bryan G. L., Norman M. L., 1998, ApJ, 495, 80 
Cavaliere A., Lapi A., Menci N., 2002, ApJL, 581, LI 
Churazov E., Sunyaev R., Forman W., Bohringer H., 2002, 

MNRAS, 332, 729 
Ciotti L., Ostriker J. P., 1997, ApJL, 487, L105 
Di Matteo T., Croft R. A. C, Springel V., Hernquist L., 

2003, ApJ, 593, 56 
Di Matteo T., Springel V., Hernquist L., 2005, Nature, 433, 

604 

Ferrarese L., Merritt D, 2000, ApJL, 539, L9 
Gebhardt K., Bender R., Bower G., Dressier A., Faber 
S. M., Filippenko A. V., Green R., Grillmair C, Ho L. C, 
Kormendy J., Lauer T. R., Magorrian J., Pinkney J., 
Richstone D, Tremaine S., 2000, ApJL, 539, L13 
Granato G. L., De Zotti G., Silva L., Bressan A., Danese 

L., 2004, ApJ, 600, 580 
Hairing N., Rix H., 2004, ApJL, 604, L89 
Harms R. J., Ford H. C, Tsvetanov Z. L, Hartig G. F., 
Dressel L. L., Kriss G. A., Bohlin R., Davidsen A. F., 
Margon B., Kochhar A. K., 1994, ApJL, 435, L35 
King A., 2003, ApJL, 596, L27 
Komatsu E., Seljak U., 2001, MNRAS, 327, 1353 
Lapi A., Cavaliere A., Menci N., 2005, ApJ, 619, 60 
Macchetto F., Marconi A., Axon D. J., Capetti A., Sparks 

W., Crane P., 1997, ApJ, 489, 579 
Marconi A., Hunt L. K, 2003, ApJL, 589, L21 
Merritt D, Ferrarese L., 2001, ApJ, 547, 140 
Murray N., Quataert E., Thompson T. A., 2005, ApJ, 618, 
569 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 
493 

Omma H., Binney J., 2004, MNRAS, 350, L13 
Quilis V., Bower R. G., Balogh M. L., 2001, MNRAS, 328, 
1091 

Reynolds C. S., Heinz S., Begelman M. C, 2001, ApJL, 
549, L179 

Reynolds C. S., Heinz S., Begelman M. C, 2002, MNRAS, 
332, 271 

Roychowdhury S., Ruszkowski M., Nath B. B., Begelman 

M. C, 2004, ApJ, 615, 681 
Ruszkowski M., Begelman M. C, 2002, ApJ, 581, 223 
Ruszkowski M., Briiggen M., Begelman M. C, 2004, ApJ, 

615, 675 

Silk J., Rees M. J., 1998, A&A, 331, LI 
Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 
Tabor G., Binney J., 1993, MNRAS, 263, 323 
Tremaine S., Gebhardt K, Bender R., Bower G., Dressier 
A., Faber S. M., Filippenko A. V., Green R., Grillmair 
C, Ho L. C, Kormendy J., Lauer T. R., Magorrian J., 
Pinkney J., Richstone D., 2002, ApJ, 574, 740 
Tucker W., David L. P., 1997, ApJ, 484, 602 
Zanni C, Murante G., Bodo G., Massaglia S., Rossi P., 
Ferrari A., 2005, A&A, 429, 399 



